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Abstract 

Consider an arbitrary local quantum field theory with a gap or an arbitrary gapless free 
theory. We consider states in such a theory, that describe two entangled particles localized 
in disjoint regions of space. We show that in such a state, to leading order, Renyi entropies 
of spatial regions, containing only one of the particles are same as their quantum mechanical 
counterparts, after subtraction of vacuum contribution. Subleading corrections depend on 
overlap of wave functions. These results suggest that Von Neumann entropy of a spatial 
region, after subtraction of vacuum contribution, can serve as a measure of entanglement of 
indistinguishable particles in pure states. 
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1 Introduction 

What makes quantum mechanics mathematically nice is linearity. On the other hand what 
makes quantum mechanics physically somewhat counterintuitive is also linearity. The best 
manifestation of this is entanglement, non triviality of which has been noted as early as 1935 [1]. 
In last few decades it has also been realized that entanglement of quantum particles can be 
put to practical use. This is at the very heart of the subject of quantum information and 
communication. Recent years have witnessed widespread interest in another somewhat different 
aspect of entanglement, namely entanglement of spatial regions in quantum field theories. This 
has been found to be relevant in diverse areas of physics such as holography [2], [3], entropic 
c theorem [4], [5], quantum phase transitions [6], [7], [8] to name a few. 

In this paper we ask a rather simple minded question, which is possibly of quite general 
interest. We consider a state in a general local quantum field theory that can be thought of as an 
entangled state of two localized particles. Then the entanglement (in the framework of quantum 
mechanics) of the emergent particles, is clearly a relevant notion in this context. On the other 
hand in the framework of field theory, the natural notion is entanglement of spatial regions 
(containing one of the localized particles). Hence one would expect some relation between 
these two notions. To our knowledge such a relation has not been investigate previously 1 . We 

1 Renyi and Von Neumann entropies of locally excited states have been studied in [9], [10] for conformal 
field theories. But there one is mainly concerned about time evolution of entanglement. In [11] mutual Renyi 
information of locally excited entangled states was considered and it was noted that this admits a quantum 
mechanical interpretation. We show this for Renyi entropies themselves and extract somewhat different physics 
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address this question in this paper and show that such a relation exists for any local quantum 
field theory. Since field theory is a more fundamental description of elementary particles than 
quantum mechanics, this also provides a finer notion of entanglement of particles * 2 . 

Our results may have interesting bearing in black hole information paradox [12], more 
precisely on Mathur’s argument [13] and AMPS argument [14]. This is because these arguments 
are based on quantum mechanical notion of entanglement, whereas the natural framework for 
the problem is quantum field theory in curved space times. To leading order, it would not 
matter, nevertheless it may be interesting to explore small corrections due to this. 

Our paper is organized as follows: in section 2, we derive our main result (2.15), which 
is relating Renyi entropies of spatial regions to those of quantum particles. In 3 we discuss 
possible usage of our results for the problem of entanglement of indistinguishable particles. 

2 Entanglement: spatial regions vs particles 

Consider a state of the following form 

| V)=N(aO i d j + bO k d l ') |0) . (2.1) 

A is a normalization constant and Oi denotes creation operator for a mode corresponding 
to the spatial wave function (/>i(x) (which is concentrated around the point aq) and internal 
index i, Oj denotes creation operator for a mode corresponding to the spatial wave function 
02 (x) (which is concentrated around the point aq) and internal index j. Similar statements 
hold for Ok, Oi. 

For a wave function or equivalently 4>(k), the corresponding state is J -^j=(fi(k)\k), 

where | k) is the one particle state of momentum k. For free theories (both gapped and gapless) 
we have explicit construction of such states in terms of creation operators. For gapped interact¬ 
ing field theories, we still have unambiguous notion of one particle momentum eigenstates \k). 
In this case we can define “creation operator” aj. by | k) =: aj/0), with |0) being the interacting 
vacuum. Since one particle states are orthogonal to the bound states as well as continnum of 
multi-particle states, we have a*, |a) = 0 for any bound/ multi-particle state | or). Thus O, O are 
simply suitable linear superpositions of these a|,-s. However the notion of one particle states 

than these papers. 

2 In spirit, this is similar to various physical quantities, e.g. magnetic moment of electron, which have a 
coarser description in quantum mechanics and a finer one in field theory. Although for entanglement the finery 
is an artefact of the state under consideration, rather than the theory. 
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is a bit ambiguous for gapless interacting theories. This is becasue being interacting one does 
not have an explicit creation operator to start with and being gapless it is not clear how to 
extract one particle states from continuous spectrum. Hence we would exclude such theories 
form our analysis. 

We choose the following normalization 3 for O -s 

<o|o?cyo) = s tj = <o|6J6j|o). ( 2 . 2 ) 

Now if the wave functions 4>i(x) and foix) are entirely supported in disjoint regions Ai and A 2 , 
such a state describes an entangled state of two particles, distinguishable by their positions. 
This may then be described, within the framework of quantum mechanics of distinguishable 
particles, as 

^ = yi^+iff 
with the spatial wave functions and 0 2 (x) serving as particle labels. 

A natural question to ask is that what is the connection between entanglement of 2.3 and 
entanglement of the region say A 1: in the state 2.1? We answer this question in this paper for 
arbitrary local quantum held theory with a gap and for arbitrary local gapless free theory. 

Bipartite entanglement of a pure state is quantified by the Von Nuemann entropy S(p) = 
— Trplogp, where p is the density matrix of any of the parties. However for our purpose, 
Renyi entropies prove more useful. n th Renyi entropy of a density matrix p is given by 

S n {p) = — log Tr p n , (2.4) 

1 — 71 

where p is the density matrix of any of the parties. A party would stand for a paricle (or a 
collection of particles) while discussing a quantum mechanical state and a region of space while 
discussing a held theory state. 

For computational reasons it is even more advantageous to look at the following quantity 

R n (p) := e (i-n)S„(p) = Tr p n . (2.5) 

First we note that we can express R„ as 

R n (p) = Tr n [p® n E^] , (2.6) 

3 Note as O-s contain only creation operators, they (anti)commute among themselves. 


(a\i)\j) + b\k)\l)) 


(2.3) 
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with 


T ?( n ) 


II' 

fc=i 


3k,lk +1 


(2.7) 


Tr n stands for trace over ' Hf n , where H\ is the Hilbert space of the subsystem under consid¬ 
eration. The indices {ik,jk} take discreet values if we are discussing a quantum mechanical 
state and continuous values if we are discussing a field theoretic state. Note that 2.6 does not 
fix E' n ' uniquely. 

As defined above, E bd is a linear operator on Hf n . One can easily extend its action on the 
full tensor product Hilbert space (Hi ® TO) 0 " by defining 


E {n \ • , . 

Jl ,0l.. .Jn, ,0n jTi ,tti ...In 


n 

k= 1 


where ik,jk are the indices of fc-tli copy of Hi and a *., 0 are the indices of 0-tli copy of "H 2 - In 
fact if pi = Tr 2 pi 2 , then 


Rn(pi) = Tr„ [/ of 2 n E <">] 


where now trace is over (Hi <E> T^) 0 " ■ From now on we will suppress the subscripts on trace 
and assume it is clear from the context. 

For us P 12 = |T)(T| (with T) given by 2.1) and subsystem that we are interested in, is 
a spatial region Q. Thus Pi = Pn = TrQ c |\I/)(\l/|, h2 c being the complementary region of fl 
The operator Eq has been studied in great details in a recent paper [11] and we will heavily 
borrow their results. If we denote local held variable at point a as q a , then E^ reads 


4 ”’ = 


nn^n-tf • - •«?}! <*(91 

j =1 aen 6ef2 


(i') 




(n') 

Qb 




X I; 


fi c ) 


( 2 , 8 ) 


where Iq c is identity operator on H® n , Hu c being the Hilbert space associated with held 
variables living in the region h2 c . Perhaps a more comprehensible way to express E^ is the 
following 

4 ”’ = / £>4 ■ ■ ■ £>« £>4 ■ ■ ■ £> 4 . ■ ■ ■ £> 4^(4 - 4 ) • • • - 4 ) 

14 ■ ■ ■ 4>(4 ■ ■ ■ 4I®I4 C • "4J(4« ■■ ■ 4J 

= I r>4...£>« £>4 e ... ...4>(4...«|®I4„ • • • 4J(4« ■ • • «J • (2.9) 


5 


Here \(/>q) denotes the state corresponding to a smooth field configuration over region fl and 
similar statement holds for \<pn c )- These are analogs of position eigenstates in quantum me¬ 
chanics. Eq has many nice properties, as explored in [11], The one particularly useful for us 
is 

i(lpl .. .^ n \E^ ] \xi ■ ■ - Xn) = ■ (2.10) 

A derivation of this property is given in A. 

Armed with these, we set out to compute R n (pn). We consider a bosonic theory to start 
with. The main result of this section, namely 2.15 would hold both for bosonic and fermionic 
fields. The corrections to this (disussed in B) would differ though. First we note 

Rn(pn) = Tr (jT® n )(T ® n |= (y 9n \Ej? } |T® n ) . 


This contains 4 n terms, each of the following form 



Here we have ignored the normalization factor N 2n . To shorten the expressions, we introduce 
some notations 


ot^) ■■■{off =--o\ 

°™)'=--°) 

o<;)...o<;)=:0 5 

6<;>...d<;>=:0,. (2.11) 

So we have 

(o|6JoJ4” ) o p 6,(o). (2.12) 

Remember, in present notation all the O, C^-s are clusters of n creation/annihilation operators. 
After some more steps, which we describe in appendix B, to set the clutter aside, this can be 
written as 

<0|4 n) |0> x 5 P2M ... 5 Puin x S jiqi ... S jnqn + other pieces . (2.13) 
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“Other pieces” are discussed in detail in appendix B. For now we just mention the relevant 
property of these pieces. They are small when the region 0 is chosen such that (f>i(x) is mostly 
supported inside 0 and 02 ( x ) is supported mostly outside 0. Thus we concentrate on the 
leading piece. When this piece is added for all the 4 n terms, one ends up with 

Af 2 "(0|4">|0) (|a| 2 ”+|l,H = (0|4”>|0) (aVM+H 2 ) 2 “«„(/5? M ), (2.14) 

where p® M is the density matrix of the first particle in the state 2.3. If O and O are supported 
in disjoint regions, which is the case considered, N ( V |a| 2 +H 2 ) _1 - Thus 

Rn (pn(tf)) ~ Rn (MO)) X R n 

or, (M*)) - S n (pn(0)) ~ S n (p? M (VO) , (2-15) 

to leading order. “Other pieces” mentioned in 2.13 and deviation of N from (|a 2 | + |6| 2 ) _1 ^ 2 
constitute corrections to this. We mention one interestig property of one of the correction 
terms (see B for detail). Among “other pieces”, second piece of B.20 contains a factor of 
\\ k jkPk- This will give non-zero contribution if i — l or j = k or i — j or k = l in 2.1. 

A little thought would convince the reader that if we consider a more general state than 
2.1, i.e. a state of the form 

|*) = JV£o«O j „O j JO) , (2.16) 

a 

2.15 would still hold. In fact our results continue to hold when instead of a single particle each 
subsystem contains a bunch of particles localized in some region. Now O, O would contain a 
bunch of creation operators, but that does not change any of the steps. 

Another interesting case is of occupation number entanglement. In this case, one would 
consider states of the form 

\4>) = J2a n O n O n \0) . (2.17) 

n 

In occupation number eigenbasis (and assuming orthogonality of the modes occupied) the 
above state would read 

M = VlM n,n) . (2.18) 

n 

We can repeat the very same steps as in B. Corresponding to each Kronecker delta in spin 
indices in B, now we would have a Kronecker delta in occupation number. Hence for present 
purpose they play the role of spin indices and 2.15 alongwith correction terms also holds for 
occupation number entanglement. 
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3 Discussion 


These results prove useful for a problem in quantum mechanics, namely the problem of en¬ 
tanglement of indistinguishable particles 4 . If wave functions of various particles are nearly 
orthogonal they can be assigned different values of some observable and the particles can ef¬ 
fectively be labelled by the value of that observable. Afterwards one can discuss entanglement 
of those particles treating them as distinguishable particles. E.g. electrons stuck to various 
lattice sites can be labelled by the site it sits on (i.e. position) and they can be entangled 
through spin. 

When wave functions start overlaping such labelling becomes ambiguous and the inherently 
indistinguishable nature of fundamental particles become important. The primary difficulty in 
discussing entanglement in such situation is identifying distinguishable subsystems, for particles 
are no more distinguishable. It is useful to take a field theoretic view of the situation, for then 
spatial regions serve as natural subsystems. When wavefunctions are localized around different 
positions, one can choose a region containing only a single particle. The subsystem then 
corresponds to that particle. This can be made quantitative if one finds a field theoretic 
quantity S, which is nearly equal to the Von Neumann entropy of the quantum mechanical 
density matrix of that particle. Even if the wave functions overlap, S remains well defined. 
Thus provided S is finite, it would be a natural candidate for entanglement of indistinguishable 
particles. The subsystem can not clearly be associated with any particle though. 

2.15 suggests S'(pq(T)) — ^(/^(O)) is the natural candidate for S, where we have used 
notations of last section and S(p) denotes the Von Neumann entropy of the density matrix p. 
This quantity has the correct limit as the wavefunction overlap approaches zero. Given that 
Von Neumann entropy quantifies entanglement in pure states [23], S(pn( x k)) indeed quantifies 
entanglement. Further subtraction of vacuum contribution renders it finite. 

2.15 also tells that S(pu(''S)) — S'(pn(0)) would in general differ from the naive quantum 
mechanical answer S, the difference being given by the “correction terms” described in B (for 
Renyi entropies). It should be noted that these corrections are never exactly zero, since there 
are no localized particles in a local quantum field theory [24], We interpret these corrections as 
field theoretic corrections to naive quantum mechanical notion of entanglement. Due to these 

4 Many attempts have been made in recent years to generalize the notion of entanglement for indistinguishable 
particles [15], [16], [17], [18], [19], [20]. Look at [21], [22] for recent review and references there in. These 
approaches mostly explore two directions, namely tensor product structure of the Hilbert space and occupation 
number entanglement between various modes (which is indeed well suited for many systems considered in 
laboratory). 



corrections, one would in general require an infinite dimensional density matrix to account for 
the vaccum subtracted Von Neumann and Renyi entropies. In the spectrum of this density 
matrix held theoretic corrections correspond to infnitely many small eigenvalues. In terms of 
the entanglement spectrum 5 , this means held theoretic effects are encoded in “high energy 
states”, i.e. quantum mechanics serves as a “low energy effective theory” in the context of 
entanglement! 

These correction terms also lead to interesting properties for the lattice analog of S, namely 
site entanglement [15], [16]. E.g. a single electron with wavefunction supported at more 
than one lattice site looks like an entangled state 6 7 ! However bizarre, such entanglement has 
actually been used in teleportation [29] and therefore physical. Given this it is natural to 
wonder whether the continuum analog of this can have some practical use as well. E.g. due 
to the held theoretic corrections, entanglement of a “single” electron can exceed log 2! It 
would be interesting to explore possible implications of this for quantum information and 
communication'. 

Acknowledgements: I am indebted to my friend Pinaki Banerjee for collaboration in earlier 
part of this work and many subsequent discussions. My special thanks to Ashoke Sen for 
various illuminating discussions and valuable comments on the manuscript. I would also like 
to thank Nilay Kundu, Arunabha Saha and whole QIC group at HRI, especially H.S. Dhar, 
M.N. Bera, Aditi Sen(De) and Avijit Misra for useful conversations and comments. Lastly I 
thank people of India for their generous support for research in theoretical physics. This work 
has been supported by DAE project 12-R&D-HRI-5.02-0303. 

5 Modular Hamiltonian H for a given density matrix p is defined as p = e~ H . Previously the spectrum of 
H 1 called “entaglement spectrum”, has been argued to carry footprints of topological order [25], in the context 
of fractional quantum Hall effect. The low energy states of the entanglement spectrum become gapless, as the 
system becomes topologically ordered. 

6 Another strange feature is that the entanglement (through spin) of two electrons at same site (distinguished 
by orbitals quantum number) is invisible to site entanglement. These points have previously been noted 
in [15], [26]. Also look at [27], [28]. 

7 We thank Arunabha Saha for pointing out this possibility. 
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A Derivation of 2.10 


To make the logic more tractable we derive 2.10 for n = 3. 

Eq ^ = J DcPqDc^qDi^qDcPq 5((fra — 4>n)S(<fra — 0n)<5(0n — 0o)l0n0n0n)(0n0n0nl 

= I Dcfr^DcfrlDfrl Dcfr^DcfrlDcfr^ I0n0o0n}(0n0n0nl x I0n>n>n c )<0n>n>nj • 

Now we can express any state \ifr) as 

10) = f DfraDfr>n c ifr[(frn,<frn c }\(fra) ® \M ■ 

Hence 

( 0 i 0 2 0 3 |£n 3) | *1X2X3) = j Dcfr^DfrlDcfrl Dcfr^Dfr^D^ 

Vbl^rn^nJ Xi[0n>0nj X2[0n>0nJ X3[0fn0nJ 

[DfoDfoDfa 

X*A<l>\iA\i c ] X* 2 [(frl,(frl c } X* 3 [(fr 3 n,(frnM(frlAh c ] V^n^nJ 03[0n>0nJ 
[ D<foD<&D<& 

Xl [0n e .0n] X3^n c >n] xt[0n c ,0n] ^l^h^n] tM0n c ,0f>] 03[0n c >0n] 
(redefining 0^ as and thinking of 0 q c as the region of interest.) 

= (X2XsXl l^l01020s)* • 

In writing ifr[frn, cfrn c ] as 0[0 a c , 0q] we are just thinking Ha ® Ha c as Ha c ® Ha, meaning now 
f2 c is our region of interest rather than fi. 

Similar derivation follows for general n and we have 

(01. • .ifrnlE^lxi ...Xn) = (X2 ■ ■ .Xi\Ea ] Mi ■ ■■i'n)* ■ (A.19) 


B Correction terms 


Now we discuss the leftover “other pieces” in 2.13. 


(a\o]o\E^o p o,\a) = <o|o'b™o p o]o,|o) + <o| 


tp(")/ 


nt nt /?("■) n 
L n i w i - C/ a 


O q \ 0 ) • 


(B.20) 
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1 st piece in B.20 can be further simplified as 


(0\0\E^O p d]0,\0) = (Q\0\Eq' > O v \a) (a|OjO g |0), (B.21) 


where we have inserted a complete set 8 of states |a)(a|. Only the \a) = |0) term survives, 
since for all others terms we have (a|O|O g |0) = 0. 

We concentrate on the piece (Q\OlE ( f^O p \Q). 



Gathering all these, we see 


(0|dj0, t 4” ) 0 r 6,|0) = (0|4”>|0) x 8 nM ... 8 nfn x S jiq , ... S iM „ 


>? ) ) t ... (o*; 1 )' 


pi 


4b Oi 


|0)‘ X « 


nqi 


jnQn 


dJ,oJ4">o p o,|o>. 


(B.22) 


In right hand side of B.22, apart from the first piece, all other pieces contain commutators 
between fields that are mostly supported in disjoint regions of space and hence are small. 
Thus the first term is the leading piece. 


C Fermionic fields 

Now we have to keep track of some signs and have to arrange stuff in terms of anti-commutators. 
We treat the cases of even and odd n separately. 

8 For interacting gapless theories there could be confusion regarding which of {|a)} represent one particle 
states. But we are excluding such theories from our analysis. 
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Even n : A generic term is still given by B.22, so all conclusions remain the same, in particular 
2.15. Only change is that end of the day we have to break up various commutators in terms 
of anti-commutators. This can easily be done using the 

n 

[A B 1 ...B n \ = ^(-l) p -'5, ... B P ^{A, B P }B P+1 ...B n . (C.23) 

p =i 


Odd n: In this case we use anti-commutators everywhere, because then we can use the 
folowing identity to break everything up into anti-commutators 

n 

{A, B 1 ... B n } = ]T(-1 y~ l B l ... Bp^iA, B P }B P+1 ...B n . (C.24) 

p =i 

We break up a generic term as 

(o|ojo;4">o r o,|o> = -{o|o;4">o r 6]o,|o) + (0|{oJ,oK>o p } O,|0). (C.25) 

In the first piece of C.25 we insert complete basis and end up with (0|OjE^O p |0)(0|OjO,j|0). 
In this, 


(0|0*4" , 0 P |0) = <o| (oW)'... o!. 1 ’ ■ ■ ■ otjw 

= —<0| (OW)* • • • (o«)' 0<;>... 0?J 4".’|0>* 

+ <o|o>”>... o«;> {4">, o‘;>... o‘:»} o) 

= -<0|4 n) |0) + <o|o<")... om {4”>, o<;>...o<:>} |0>. 

(C.26) 


Gathering all these, we see 


<o|6Jo, t 4">o p 6,|o> = (o|4">|o> s n „...s inPI x s jm ...s jntL 

- <oio<")... ow o<:>... o<:»} io)^ ... s jn „ 

+ (0|{o],0. t 4” l 0 p }o,|0) (C.27) 


Since the leading piece remains the same, 2.15 goes through. The details of the corrections 
terms is different from those in B, although the difference is only quantitative. 
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